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A  CFD  study  was  conducted  for  the  electrolyte  transport  in  a  single-cell  alkaline  fuel  cell.  The  transport 
model  that  accounts  for  the  balance  of  mass,  momentum,  energy,  species  and  electric  charge  in  the 
electrolyte  was  solved.  The  CFD  results  captured  the  presence  of  primary  and  secondary  (shunt)  currents 
in  the  liquid  electrolyte.  The  shunt  currents  were  present  in  the  regions  adjacent  to  the  separator  inlet 
and  exit,  but  not  inside  the  separator  nor  in  the  extended  channels.  The  calculated  local  shunt  current 
density  in  the  corner  regions  of  the  separator  inlet  and  exit  was  higher  than  the  average  primary  current. 
The  results  also  showed  that  the  IR  loss  due  to  the  electric  potential  drop  at  high  current  densities  is 
comparable  to  the  overpotential  loss  at  the  anode.  Passive  and  active  mitigation  measures  against  the 
shunt  currents  were  suggested. 
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1.  Introduction 

A  stationary  fuel  cell  is  a  viable  option  for  electrical  power  gen¬ 
eration  from  renewable  sources.  For  example,  H2  produced  by  solar 
or  wind  turbine  based  electrolysis  can  be  used  to  fuel  the  fuel 
cells.  Among  the  fuel  cell  platforms,  the  alkaline  fuel  cell  (AFC) 
remains  a  good  prospect  as  it  has  high  reliability  and  it  can  use 
non-precious  metal  for  the  electrodes  [1-7].  Flowever,  to  be  com¬ 
petitive  with  other  power  sources,  the  power  density  (as  well  as 
the  cost  of  manufacturing,  operation  and  maintenance)  needs  to 
be  improved.  In  [8],  it  is  shown  that  maintaining  a  uniform  elec¬ 
trolyte  concentration  is  an  effective  strategy  to  increase  the  limiting 
current  density.  Another  technology  challenge  is  corrosion  due  to 
the  presence  of  shunt  currents.  Mathematical  models  describing 
the  electrochemical  and  transport  processes  in  AFC  electrolytes  are 
available;  for  example,  see  [9-15].  However,  the  transport  models 
were  formulated  assuming  ideal  solutions.  In  [8]  it  is  shown  that 
the  ideal  solution  formulation  underestimates  the  limiting  current 
density  by  50%  for  the  concentrated  electrolyte  condition  consid¬ 
ered,  7M.  Furthermore,  the  effects  of  temperature  on  transport 
properties  were  also  neglected,  which  was  shown  to  underesti¬ 
mate  the  limiting  current  density  by  30%  when  compared  to  an 
adiabatic  boundary  condition  [8].  In  this  paper,  the  ideal  solution 
and  isothermal  assumptions  were  relaxed  and  a  CFD  calculation 
was  performed  for  prediction  of  the  presence  of  shunt  currents  in 
alkaline  fuel  cells. 
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Shunt  currents  or  bypass  currents  are  present  in  liquid  elec¬ 
trolytes,  coolant  circuits  of  fuel  cell  stacks,  and  bipolar  electrolyzers 
[6,1 6-23  ].  The  presence  of  shunt  currents  is  the  result  of  an  electric 
potential  gradient  along  a  common  electrolyte  path  that  forms  a 
driving  force  for  the  net  ionic  fluxes.  Shunt  currents  are  not  desired 
as  they  lead  to  corrosion  problems,  reduced  service  life,  and  sys¬ 
tem  power  losses.  Shunt  currents  can  be  reduced  or  eliminated  by 
applying  (a)  a  protective  current  against  the  shunt  current  direc¬ 
tion,  (b)  protective  coatings  on  surfaces  of  the  electrolyte/coolant 
channels,  or  (c)  baffle  plates  to  stop  the  ionic  fluxes.  Understanding 
the  transport  phenomena  in  electrolytes  can  help  design  consid¬ 
erations  for  effective  mitigation  measures.  To  model  the  current 
transport,  models  based  on  equivalent  electrical  circuits  [24-28] 
or  Laplace  equations  [29-32]  are  developed.  To  further  the  state 
of  knowledge,  a  transport  model  that  accounts  for  the  balance  of 
mass,  momentum,  energy,  species,  and  charge  in  electrolytes  and 
electrodes,  and  that  satisfies  the  thermodynamic  constraints  of  H2 
and  02  dissolving  in  electrolytes  is  implemented  in  this  study.  The 
objective  of  this  paper  is  to  conduct  a  CFD  based  study  of  the  elec¬ 
trochemical  and  transport  processes  in  electrolytes,  to  gain  insights 
into  the  presence  of  shunt  currents  in  AFCs,  and  to  provide  sugges¬ 
tions  that  aid  in  design  considerations. 

2.  Formulation 

The  assumptions  invoked  in  the  formulation  are  (a)  Newto¬ 
nian  fluid,  (b)  negligible  pressure  effects  on  enthalpy,  (c)  negligible 
magnetic  effects  due  to  electric  field,  (d)  negligible  pressure  and 
temperature  effects  on  diffusion  transport,  (e)  negligible  vis¬ 
cous  dissipation,  pressure  work,  and  Dufour  effects  in  energy 
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Nomenclature 

C  Molar  concentration  (kmol  m-3 ) 

F  Faraday  constant  (C  kmol-1 ) 

i  Current  density  (A  m-2 ) 

i  Vector  quantity  of  current  density  (A  m-2 ) 

Ru  Universal  gas  constant  (J  kmol-1  K-1 ) 

Slpa  Source  term  of  Eq.  ( 1 ) 

T  Temperature  (I<) 

t  Transference  number 

x  Coordinate,  cross-stream  direction  (Fig.  1 ) 

y  Coordinate,  streamwise  direction  (Fig.  1 ) 

Greek 

8  Dimensionless  distance;  5  =  0  and  1  denotes  separa¬ 

tor  inlet  and  exit,  respectively 
A0  Electric  potential  drop  (V) 

rj  Overpotential  (V) 

Y±,c  Mean  activity  coefficient  based  on  molarity 

k  Electrical  conductivity  (S  m-1 ) 

kd  Diffusional  ionic  conductivity  (A  m-1 ) 

0  Electrical  potential  at  electrolyte  phase  (V) 

Subscripts 
a  Anode 

c  Cathode 

e  Electrolyte 

j5  Liquid  phase 

o  Solid  phase 

0  Solvent 

-  Anion 

Superscripts 
eff  Effective 

D  Diffusion 

i  Charge 

0  Reference  state 


equation,  (f)  no  homogeneous  chemical  reactions,  (g)  electroneu¬ 
trality,  (h)  homogeneous  and  continuous  media  for  gas  and  liquid 
phase,  (i)  local  thermal  equilibrium,  (j)  vapor-liquid  equilibrium 
at  the  gas-liquid  interface,  and  (k)  negligible  electrical  resistance 
of  the  electrodes.  The  governing  equations  describing  the  trans¬ 
port  of  mass,  momentum,  species,  energy  and  charge  balance  are 
given  in  [8],  which  account  for  (a)  mass  and  species  addition  or 
removal  due  to  gas  solubility  and  H20  phase  change,  (b)  momen¬ 
tum  exchange  due  to  Darcian  flow  in  porous  electrodes,  (c)  energy 
source  or  sink  terms  due  to  heat  and  entropy  generation  at  cata¬ 
lyst  layers,  Joule  heating  in  separator  and  catalyst  layers,  and  latent 
heat  due  to  H20  phase  change,  and  (d)  charge  generation  at  cata¬ 
lyst  layers.  The  constitutive  equations  describing  the  reaction  rate 
(Butler-Volmer  equation),  capillary  pressure  (Leverett-J  function), 
and  physical  properties  are  given  in  [8].  Species  equation  of  the  liq¬ 
uid  phase  is  solved  only  for  OH-.  The  K+  concentration  is  calculated 
from  electroneutrality. 


i 


1:  Anode  Co*  Channel 
2:  Anode  Cai  Diffusion  Layer 
3:  Anode  Caiafytt  La^er 
4:  Separator  Lajer 
5:  Cathode  Catalyst  Layer 
6:  Cathode  Cas  Diffusion La^ex 
7:  Cathode  Cat  Channel 


KOH 


Fig.  1.  Schematic  of  alkaline  fuel  cell. 


software  is  utilized.  The  user  defined  functions  (UDFs)  were  used 
for  calculation  of  the  physical  properties  and  the  source  terms  of  the 
governing  equations,  and  to  satisfy  the  prescribed  boundary  con¬ 
ditions.  The  graphical  user  interface  (GUI)  was  used  to  input  the 
UDFs  to  the  corresponding  governing  equations.  The  user  defined 
scalar  (UDS)  was  introduced  to  solve  the  charge  equation.  Two  con¬ 
vergence  criteria  were  enforced  to  ensure  that  (a)  the  local  cell 
voltage  is  the  same  everywhere  along  the  separator,  and  (b)  the 
total  current  density  equals  that  calculated  by  the  specified,  aver¬ 
aged  cell  current  density.  Description  of  the  computational  method 
and  validation  of  the  CFD  results  are  given  in  [8]. 

4.  Results  and  discussion 

The  primary  and  secondary  (shunt)  currents  are  calculated  by 
charge  balance  equation  [8]: 

V  .  V<fy)  +  V  .  (^effV  In  Cpe)  =  (1 ) 

The  source  term  in  the  charge  equation  is  used  to  calculate  the 
charge  transfer  between  the  solid  and  liquid  phase.  The  current 
density,  i,  is  calculated  from  [8]: 

i=  -kV0-  kd  VlnCe  (2) 

where  k  is  the  electrical  conductivity,  0  the  electric  potential,  Ce 
the  electrolyte  molar  concentration,  and  kd  the  diffusional  ionic 
conductivity  [33]  or  diffusion  conductivity  in  short.  The  diffusion 
conductivity,  which  is  associated  with  the  ionic  transport  due  to 
the  presence  of  a  concentration  gradient,  is  calculated  by: 


3.  Solution  method 

The  computational  domain  was  divided  into  seven  regions 
(domains)  as  illustrated  in  Fig.  1.  Domains  1,  2,  and  3  are  the 
anode  gas  channel,  gas  diffusion  layer  (GDL),  and  catalyst  layer, 
respectively.  Domain  4  is  the  separator.  Domains  5,  6,  and  7  are 
the  cathode  catalyst  layer,  GDL,  and  gas  channel,  respectively. 
To  solve  the  coupled  nonlinear  governing  equations,  FLUENT® 


(2RuTk\  f  l-t°+Ce\  / 1  +d\ny±c  \ 
1  F  i\  2 C0  J  \  dlnCe  J 


(3) 


where  Ru  is  the  universal  gas  constant,  T  the  temperature,  F  the 
Faraday  constant,  t°_  the  reference  state  transference  number  of 
negative  ion,  and  y±<c  the  mean  molar  activity  coefficient.  Con¬ 
stant  and  uniform  electric  potential  is  assumed  for  the  electrodes; 
namely  0a  =  O  and  Ecen  for  the  anode  and  cathode,  respectively. 
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(b) 

Fig.  2.  KOH  concentration  along  the  electrodes  (5  =  0  and  1  denotes  separator  inlet 
and  exit,  respectively)  at  specified  average  current  density,  adiabatic  boundary  con¬ 
dition;  (a)  anode  and  (b)  cathode. 

Electrolyte  concentration  and  electric  potential  are  solved  in  the 
CFD  calculation.  The  primary  and  shunt  currents  are  calculated 
accordingly. 

For  the  computation,  the  baseline  condition  given  in  Table 
6  of  [8]  was  used  unless  noted  otherwise.  In  short,  the  KOH 
concentration  and  temperature  were  set  to  7  M  and  80  °C  at  the  sep¬ 
arator  inlet,  respectively.  The  system  pressure  was  set  to  4.1  atm, 
and  with  pure  H2  and  02  as  the  fuel  and  oxidant.  The  relative 
humidity  of  H2  and  02  was  set  to  0%.  The  anode,  cathode  and 
electrolyte  velocities  at  the  separator  inlet  were  set  to  0.2,  0.1, 
and  0.01  ms-1,  respectively.  An  adiabatic  boundary  condition  was 
imposed. 

41.  Electrolyte  concentration 

The  electrolyte  concentration  distributions  along  anode  and 
cathode  catalyst  layers  at  specified  current  densities  are  shown  in 
Fig.  2(a)  and  (b),  respectively.  At  fixed  streamwise  locations  (5),  the 
KOH  concentration  at  anode  decreases  with  increasing  current  den¬ 
sity.  This  trend  is  expected  as  OH-  is  consumed  at  anode.  The  KOH 
concentration  was  below  the  specified  7  M  at  the  separator  inlet 


0.00  0.20  0.40  0.60  0.80  1.00 

Dimensionless  Streamwise  Distance 

Fig.  3.  Local  current  density  along  the  electrodes  at  specified  current  density 
i  =  1.0  A  cm-2,  adiabatic  boundary  condition. 


for  all  current  densities  examined  except  for  low  current  densities, 
e.g.,  i<  0.2  A  cm-2.  For  a  fixed  current  density,  the  KOH  concentra¬ 
tion  decreases  along  the  streamwise  direction  for  i<  1.0  A  cm-2. 
However,  for  high  current  densities,  e.g.,  i>  1.0  A  cm-2,  the  KOH 
concentration  decreases  sharply  immediately  downstream  of  the 
separator  inlet  and  gradually  increases  toward  the  separator  exit, 
consistent  with  the  trend  of  the  calculated  local  current  density 
shown  in  Fig.  3.  For  the  cathode,  where  OH-  is  produced,  the 
KOH  concentration  at  fixed  streamwise  location  increases  with 
increasing  current  density.  Along  the  cathode,  the  predicted  KOH 
concentration  was  above  the  specified  electrolyte  concentration 
7  M  at  the  separator  inlet.  The  results  showed  that  it  first  increased 
and  then  decreased  gradually  along  the  cathode.  The  increase  is 
due  to  the  combined  effects  of  high  local  current  density  at  sep¬ 
arator  inlet,  e.g.,  see  Fig.  3,  and  loss  of  water  to  dry  02  stream  in 
the  cathode  gas  channel.  For  i  =  0.0  A  cm-2,  loss  of  H20  to  dry  gas 
channel  due  to  evaporation  results  in  KOH  concentrations  higher 
than  7  M.  The  results  suggest  that  for  low  current  densities  the 
gas  phase  transport  due  to  the  loss  of  H20  to  gas  streams  is  the 
predominant  mechanism  that  determines  the  KOH  concentration 
along  the  anode  and  cathode,  whereas  for  high  current  density  con¬ 
ditions,  the  reactant  dissolving  rates,  the  transport  mechanisms 
and  the  electrochemical  reaction  rates  are  the  predominant  mech¬ 
anisms. 

4.2.  Electric  potential 

The  electric  potential  drop  across  the  separator  (liquid  elec¬ 
trolyte)  at  specified  current  densities  is  shown  in  Fig.  4.  The  results 
show  that  the  electric  potential  drop  increases  with  increasing  cur¬ 
rent  density,  and  it  reaches  0.07  V  at  the  limiting  current  density 
2.3  A  cm-2.  In  fact  at  2.3  A  cm-2  the  electric  potential  drop  across 
the  electrolyte  is  comparable  to  that  of  the  overpotential  loss  at 
the  anode.  More  discussion  on  overpotential  losses  will  follow. 
The  simulation  also  shows  that  for  high  current  densities  (e.g., 
i>l  A  cm-2)  the  electric  potential  drop  jumps  to  a  higher  value 
immediate  downstream  of  the  separator  inlet.  The  potential  drop 
then  decreases  slightly  at  downstream  locations.  For  low  current 
densities,  the  electric  potential  drop  decreases  slightly  along  the 
streamwise  direction. 
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Fig.  4.  Liquid-phase  electric  potential  drop  in  KOH  solution  as  functions  of  current 
density  (i)  and  streamwise  location  (5),  adiabatic  boundary  condition. 


The  distribution  of  electric  field  potential  in  the  liquid  elec¬ 
trolyte  is  shown  by  the  contour  plot  given  in  Fig.  5;  the  current 
density  is  set  to  1 A  cm-2.  It  is  noted  that  the  presence  of  the  elec¬ 
tric  field  potential  gradient  in  the  streamwise  direction  leads  to 
the  development  of  shunt  currents.  The  local  primary  current  den¬ 
sity  decreases  along  the  anode  and  cathode,  but  it  remains  the 
same  at  the  anode  and  cathode  for  fixed  streamwise  locations 
(cf.  Fig.  3).  The  decrease  of  current  density  along  the  streamwise 
direction  is  consistent  with  the  decrease  of  reactant  concentration 
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Fig.  5.  Contour  plot  of  liquid-phase  electric  potential  in  KOH  solution  at  specified 
current  density,  i  =  1 A  cm-2,  adiabatic  boundary  condition.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
the  article.) 


(a) 


(b) 

Fig.  6.  Local  overpotential  at  electrodes  as  functions  of  current  density  (i)  and 
streamwise  location  (5),  adiabatic  boundary  condition;  (a)  anode  and  (b)  cathode. 


Fig.  7.  Net  shunt  current  density  as  a  function  of  streamwise  direction  (5)  at  speci¬ 
fied  current  density  i  =  1 A  cm-2,  adiabatic  boundary  condition. 
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Fig.  8.  Contour  plot  of  shunt  current  density  inside  the  separator  (0.05  <y/m  <  0.10)  and  extended  channel  at  specified  current  density  i  =  l  Acm-2,  adiabatic  boundary 
condition;  green  color  indicating  zero  shunt  current.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


(i.e.,  H2  and  02)  along  the  gas  channels,  suggesting  that  the  species 
(i.e.,  OH-)  transport  across  the  electrolyte  is  the  predominant 
mechanism  responsible  for  setting  the  local  primary  currents. 

Distributions  of  the  local  overpotential  at  anode  and  cathode 
are  shown  in  Fig.  6(a)  and  (b).  The  anode  overpotential  increases 


with  increasing  current  density.  But  the  overpotential  remains  rela¬ 
tively  uniform  in  the  streamwise  direction.  At  low  current  densities 
(e.g.,  i<lAcm-2),  the  overpotential  is  highest  at  the  separator 
inlet  (i.e.,  5  =  0),  which  is  a  result  of  the  high  local  current  den¬ 
sity  predicted  (Fig.  3 ).  At  high  current  densities  (e.g.,  i  =  2.3  A  cm-2 ), 


(a)  Separator  Outlet  Adiabatic  Boundary  (b)  Separator  Inlet 


Fig.  9.  Vector  plot  of  current  density  inside  separator  and  extended  channel  at  specified  current  density  i  =  1 A  cm-2,  adiabatic  boundary  condition;  vector  length  representing 
magnitude  of  the  shunt  current  density. 
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(a) 


Streamwise  Distance  /  m 

(b) 


Streamwise  Distance  /  m 

(C) 

Fig.  10.  Net  shunt  current  density  versus  streamwise  distance  inside  the  separator  (0.05  and  0.10m)  and  in  extended  channel  at  specified  current  density  i  =  l  Acrrr2; 
isothermal  boundary  condition:  (a)  80  °C,  (b)  50  °C,  and  (c)  25  °C. 


the  overpotential  first  decreases  then  increases  along  the  stream- 
wise  direction,  see  Fig.  6(a).  Similar  trends  were  predicted  for 
the  cathode  overpotential,  except  for  high  current  densities  (e.g., 
i  =  2.3  A  cm-2 ).  At  such  conditions,  the  overpotential  first  increases 
then  deceases  toward  the  separator  exit.  It  is  noted  that  the  cathode 
overpotential  was  significantly  higher  than  the  anode  overpoten¬ 
tial,  0.4  V  versus  0.08  V,  which  is  consistent  with  the  well  known 
fact  that  the  cathode  overpotential  is  the  predominant  loss  mech¬ 
anism.  The  low  cathode  exchange  current  density  is  responsible 
for  the  high  overpotential  at  low  current  densities.  At  high  current 
densities,  availability  of  the  reactants  at  cathode  sets  the  overpo¬ 
tential  losses,  i.e.,  it  is  dominated  by  the  cathode  concentration 
overpotential. 

4.3.  Shunt  current 

4.3 A.  Local  and  net  shunt  current 

The  net  shunt  current  across  the  separator  that  accounts  for  the 
net  ion  fluxes  in  the  streamwise  direction  is  shown  in  Fig.  7.  As 
shown,  very  low  net  shunt  currents  are  predicted  inside  the  sepa¬ 
rator  and  in  the  extended  channels.  Two  peaks  are  predicted  at  the 
separator  inlet  and  exit,  i.e.,  at  streamwise  distances  of  0.05  and 


0.10  m,  respectively.  The  calculated  net  shunt  current  densities  are 
0.09  and  0.11  Acm-2  at  the  inlet  and  exit.  The  color  coded  local 
shunt  current  contours  are  plotted  in  Fig.  8.  The  results  showed  that 
the  local  shunt  current  density  in  the  corner  regions  of  the  separator 
inlet  and  exit  was  as  high  as  1.8  Acm-2.  These  large  shunt  current 
densities  are  to  satisfy  the  electric  potential  equation  in  the  regions 
adjacent  to  the  separator  inlet  and  exit.  The  vector  plot  of  the  cur¬ 
rent  density  is  shown  in  Fig.  9.  The  expanded  views  at  the  separator 
exit  and  inlet  are  shown  by  inserts  (a)  and  (b)  in  Fig.  9.  At  the  sepa¬ 
rator  inlet  and  exit,  shunt  currents  are  seen  leaving  the  anode  and 
returning  to  cathode.  No  current  is  observed  in  the  extended  chan¬ 
nel.  In  the  regions  where  shunt  currents  are  present,  electrolysis 
cells  may  form  based  on  the  local  electric  potential,  pH  value  and 
material  properties  that  lead  to  corrosion  problems  [34].  The  results 
suggest  that  mitigation  measures  are  needed  to  address  potential 
corrosion  problems  near  the  separator  inlet  and  exit. 

To  examine  the  effects  of  thermal  boundary  conditions  on  the 
presence  of  shunt  currents,  three  isothermal  boundary  conditions 
were  examined:  (a)  80  °C,  (b)  50  °C,  and  (c)  25  °C.  As  shown  in 
Fig.  10,  the  peak  shunt  current  density  varies  with  varying  thermal 
boundary  conditions.  The  local  shunt  current  density  distribu¬ 
tion,  however,  remains  qualitatively  similar  for  the  adiabatic  and 
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Fig.  11.  Contour  plot  of  shunt  current  density  of  isothermal  boundary  conditions  (25, 50,  or  80  °C)  at  specified  current  density  i  =  1 A  cm-2 ;  (a)  separator  inlet  and  (b)  separator 
outlet.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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isothermal  boundary  conditions  considered,  cf.  Figs.  8  and  11.  This 
finding  suggests  that  shunt  currents  will  be  present  regardless  of 
the  thermal  boundary  condition,  and  that  thermal  management 
will  not  be  an  effective  means  to  mitigate  the  formation  of  shunt 
currents  although  it  has  a  significant  effect  on  the  electrolyte  trans¬ 
port  [8]. 

4.3.2.  Shunt  current  mitigation 

The  presence  of  large  electric  potential  gradients  in  the  separa¬ 
tor  inlet  and  exit  regions  is  responsible  for  the  presence  of  shunt 
currents,  and  potentially  for  forming  electrolysis  cells  around  the 
electrodes.  The  electrochemical  redox  reaction  in  the  electroly¬ 
sis  cells  is  a  source  of  corrosion  [28,34].  To  reduce  the  corrosion, 
passive  and  active  preventive  methods  are  developed.  The  pas¬ 
sive  method  includes  the  application  of  protective  coatings  to  the 
electrode  surface  in  contact  with  the  electrolyte,  or  placement  of  a 
porous  barrier  in  electrolyte  passages  to  increase  the  resistance  in 
the  regions  of  concern  or  to  capture  the  net  ionic  fluxes  leaving  the 
separator.  A  review  of  the  coating  materials  can  be  found  in  liter¬ 
ature  [35],  and  porous  barriers  or  shields  to  prevent  shunt  current 
corrosion  are  reported  in  US  patents  [36-39].  The  simulation  results 
suggest  that  passive  measures  be  applied  in  the  separator  inlet 
and  exit  regions.  For  the  active  protective  methods,  the  simulation 
suggests  an  opposing  voltage  that  counters  the  electric  potential 
gradient  in  the  electrolytes  such  as  those  reported  in  [28,40,41] 
should  be  applied  immediately  upstream  of  the  separator  inlet  and 
downstream  of  the  separator  exit. 

5.  Summary  and  conclusion 

A  CFD  study  was  conducted  for  the  electrolyte  transport  in  a 
single-cell  alkaline  fuel  cell.  The  transport  model  accounts  for  the 
balance  of  mass,  momentum,  energy,  species  and  electric  charge 
in  the  electrolyte,  and  was  solved  using  a  commercial  CFD  soft¬ 
ware  FLUENT®.  The  user  defined  functions  were  developed  for 
calculation  of  the  physical  properties  and  the  source  terms  of  the 
governing  equations,  and  to  satisfy  the  prescribed  boundary  con¬ 
ditions. 

The  CFD  study  captured  the  presence  of  the  primary  and  sec¬ 
ondary  (shunt)  currents  in  the  electrolyte  of  a  single-cell  alkaline 
fuel  cell.  The  primary  currents  varied  along  the  anode  and  the 
cathode  and  were  highest  at  the  separator  inlet.  The  maximum 
local  current  density  was  more  than  50%  higher  than  the  specified 
average  current  density.  And  identical  local  current  densities  were 
calculated  for  the  anode  and  cathode  everywhere  but  the  separator 
exit  where  shunt  currents  were  present.  The  following  conclusions 
can  be  drawn: 

1.  The  predicted  local  shunt  current  is  in  the  regions  where  there 
exists  a  large  electric  potential  gradient  and  where  the  primary 
current  is  lacking,  i.e.,  in  the  regions  adjacent  to  the  separator 
inlet  and  exit.  Shunt  currents  are  not  predicted  in  the  extended 
channels. 

2.  Shunt  currents  are  present  regardless  of  applied  thermal  bound¬ 
ary  conditions  being  adiabatic  or  isothermal  (25,  50,  and  80  °C). 

3.  The  level  of  humidification  of  the  reactant  gas  streams  is  impor¬ 
tant  to  KOH  concentration  at  low  current  density  conditions, 
through  the  loss  of  H20  due  to  evaporation  to  gas  channels. 

4.  At  high  current  densities,  the  IR  loss  due  to  the  electric  potential 
drop  is  comparable  to  the  overpotential  loss  at  the  anode;  both 


are  about  an  order  of  magnitude  lower  than  the  overpotential 
loss  at  the  cathode. 
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